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ABSTRACT 

During the epoch of reionization the 21cm signal is sensitive to the scattering rate of the ultraviolet photons, 
redshifting across the Lya resonance. Here we calculate the photon scattering rate profile for a single ultraviolet 
source. After taking into account previously neglected natural broadening of the resonance line, we find that 
photons approach the resonance frequency and experience most scatterings at a significantly smaller distance 
from the source than naively expected r = (Av/vo)(c/H), where Av = v- is the initial frequency offset, and 
the discrepancy increases as the initial frequency offset decreases. As a consequence, the scattering rate P a (r) 
drops much faster with increasing distance than the previously assumed 1 jr 1 profile. Near the source, (r < 1 
comoving Mpc), the scattering rate of photons that redshift into the Lya resonance converges to P a (r) oc r" 7 / 3 . 
The scattering rate of Lya photons produced by splitting of photons that redshift into a higher resonance (Ly7, 
Ly<5, etc.) is only weakly affected by the radiative transfer, while the sum of scattering rates of Lya photons 
produced from all higher resonances also converges to P a (r) oc r" 7 / 3 near the source. At 15 < z < 35, on scales 
of ~0.01-20/r'Mpc (comoving), the total scattering rate of Lya photons from all Lyman resonances is found 
to be higher by a factor of ~ 1 + 0.3 [(1 +z)/20] 2 / 3 than obtained without full radiative transfer. Consequently, 
during the early stage of reionization, the differential brightness of 21cm signal against the cosmic microwave 
background is also boosted by a similar factor. 

Subject headings: cosmic microwave background - cosmology: theory - diffuse radiation - intergalactic 
medium - radio lines: general 



1. INTRODUCTION 

The next generation of radio telescopes (e.g., LOFAR, 
MWA, SKA, and 21CMA) 4 promises to open a first obser- 
vational window into the epoch preceding the end of reion- 
ization at z > 6. By measuring the redshifted 21cm signal 
from neutral hydrogen, the new telescopes can provide us 
with information on the history of reionization, the nature of 
the first radiation source, the spectrum of the primordial den- 
sity perturbation field, the cosmological paramete rs, the phys- 
ical p r operties of dark m a tter pa rt icles, etc. (e.g.jMadau et al.1 
| 1997t Iciardi & Madaul l2003t lLoeb & Zaldarriaeal 1 20041: 
Nusserl 120051: iMcOuinn et all 120061; iBarkana & Loebl 120051 
Ichuzhov et alJl2006tlShchekinov & Vasilievll2006l) . 

The 21cm signal can be observed either in absorption or 
emission against the cosmic microwave background (CMB), 
when the hydrogen spin temperature, T s , is different from the 
CMB temperature, 7cmb ■ The former is defined by the relative 
populations of the hyperfine states of hydrogen atoms 
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where = 1 .4 GHz is the frequency of hydrogen hyperfine 
transition. In the high-density gas clouds, which so far have 
been the only detected sources of the 21cm signal, the decou- 
pling of T s from 7cmb is done by collisions between atoms, 
which induce direct transitions between the hyperfine states 
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and couple the spin temperature to the hydrogen kinetic tem- 
perature, 7fc. In the intergalactic medium (IGM), on the other 
hand, after z ~ 30 the density becomes too low for collisional 
coupling to be effective and the decoupling of 7^ from 7cmb 
can be effectively done only b y scatterings of Ly a photons, 
which likewise couple T s to 7* dWoufhuvsenll 19521) . The spin 
temperature is therefore a weighted function of 7cmb and 7^ 
(lFieldlll958l) 
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where y c is the collisional coupling constant (which we ne- 
glect throughout the paper). The radiative coupling constant, 

y a , is 
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where fa = 0.4162 is the oscillator strength of the Lya tran- 
sition, r„ = hv^Jk = 0.068 K, A 10 is the spontaneous emis- 
sion coefficient of the hyperfine transition and 7o is the in- 
tensity at Lya resonance, when the backreaction on the in- 
cident photons caused by resonant scattering is neglected. 
For the unperturbed IGM the backreaction correction, S a , is 
dChuzhov & Sh apiro 20061) 
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The Lya photons can be produced in several ways, in- 
cluding recombinations, line cooling, and collisional excita- 
tion of atoms by non-thermal electrons produced by X-rays 
dChuzhov et al.ll2006t IChen & Miralda-Escudell2006l) . How- 
ever, in case the reionization epoch was dominated by stellar 
ultraviolet (UV) sources (which at present is the most popu- 
lar theory) most Lya photons in the neutral IGM originate as 
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photons between Lya and Ly-limit frequencies. Due to ex- 
pansion of the Universe, a fraction of these photons (those 
emitted between Lya and Ly/3) would gradually redshift into 
the Lya resonance. The rest would redshift into one of the 
higher resonances (generally the one which is just below their 
initial frequency). The scatterings of high resonance photons 
produce electron excitations to the np (n > 2) states, which 
are shortly followed by deexcitations and emission either of 
the original photon (in case the cascade goes directly to the 
ground state, np — > Is) or of several lower energy photons 
(in case the cascade goes via some intermediate state). Since 
between ~ 12% and 23% of cascades follow the latter route, 
after just a few scatterings most of the high resonance photons 
will be split. Depending on whether the last cascade goes via 
the 2p or 2s the resulting photons will or will not include a 
Lya photon. 

In this paper, we calculate the Lya scattering rate P a {r) as a 
function of distance r from the UV radiation source. Contrary 
to the naive expectation, the scattering rate does not evolve 
as 1/V 2 . While P a (r) oc 1 /r 2 is a good approximation at low 
redshifts when hydrogen is mostly ionized, at high redshifts 
where the Gunn-Peterson optical depth to the Lya photons 
is extremely large, we find that the scattering rate profile be- 
comes much steeper. Therefore, until the radiation intensity 
reaches the saturation levels (i.e., y a ^> 7cmb/^1)> the fluctu- 
ations of the 21cm signal would be significantly stronger than 
previously estimated. 

The paper is organized as following. In § 2, we describe 
our calculation of radiative transfer for "continuum" Lya pho- 
tons (i.e., photons that gradually redshift into the Lya res- 
onance). In § 3, we describe the radiative transfer of "in- 
jected" Lya photons (i.e., photons produced by splitting of 
high resonance photons). In § 4, we study the total scat- 
tering rate from "continuum" and "injected" Lya photons. 
In § 5, we summarize and discuss our results. Through- 
out the paper, we adopt a spatially-flat ACDM cosmologi- 
cal model with matter density il m = 0.25 and baryon density 
fib = 0.044 in units of the critical density and a Hubble con- 
stant Hq = 100k = 72kms _I Mpc _1 , which is consisten t with 
the constraints from the WMAP data (ISpergel et al.ll2007h . 

2. RADIATIVE TRANSFER - "CONTINUUM" PHOTONS 

The fate of the UV photon depends on the frequency at 
which it was emitted. Most photons originally emitted be- 
tween Lya and Ly/3 frequencies will travel a large distance, 
up to (i/0/v a -l)c/H(z) ~ 300[(1 +z)/26]" 1 / 2 Mpc (comov- 
ing), before being scattered by one of the hydrogen atoms. As 
photon redshifts closer to the Lya resonance its scatterings 
become more frequent and the mean distance between subse- 
quent scatterings rapidly drops (see Figure [TJ. Therefore al- 
most all scatterings occur within a very small region 5 . Thus to 
make an accurate estimate of the scattering rate, it is in prac- 
tice sufficient to count only scatterings occurring within this 
region, where photons are very close to the resonance and the 
number of scatterings is of the order of Gunn-Peterson opti- 
cal depth. However, to derive the scattering rate profile, P a (r), 
one also needs to calculate the distance of this region from the 
radiation source, which by contrast is determined mainly by 
the first few scatterings before photons redshift into the Lya 
resonance. 

In general, photon scattering cross-section is the function 

5 We find that at z ~ 25 above 99% of the scatterings occur within ~ 1 .5 
kpc comoving. 



of both its frequency and the gas temperature. Neglecting hy- 
perfine splitting, the cross-section without thermal broadening 
can be written as 
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with the line profile being a Lorentz function 
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where vq is the line-center frequency, / is the oscillator 
strength and 7 is the spontaneous decay rate (7/471 ~ 10~ 8 ^o 
for Lya ). Thermal broadening introduces a core in the cross- 
section around the line center with a width of the order of 
thermal velocity in units of light speed c, which is about 
10~ 6 for the unheated IGM in the redshift range we are in- 
terested. In most calculations of the Lya scattering profile 
around first sources, the line profile <j)(y) is simplified to a 
Dirac 80 function and scatterings occur as a photon redshifts 
into line-center frequency. We show in this paper that the fre- 
quency dependence of the cross-section affects the scattering 
rate profile. 

Since in the high optical depth regime, the first scatterings 
happen when photon frequency is still significantly above res- 
onance (i.e., in the Lorentz wing), the distance between the 
source and the point where photon enters into the resonance 
is nearly independent of gas temperature. Hence the normal- 
ized scattering rate profile is also independent of temperature 
[though not the total number of scatterings, see eq. Q]. Since 
the line frequency offset Av = v-vq relevant for our study 
is much larger than the quantum width and the thermal core 
width, the cross-section [eq. ©], being in the wing regime, 
can be well approximated as 
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where a c oc /7/Vq is formally the scattering cross-section at 
twice the line-center frequency. We note that in this study, 
the initial frequency of a photon is bluewar d of line center, 
i.e., A v > 0. This is different from that in Loe b & Rvbickil 
( 1999). They investigate the brightness and spectral distribu- 
tions of escaping Lya radiation around sources before reion- 
ization, and in their case photons start at a frequency slightly 
redward of the resonance, i.e., Av < 0. 

As a photon travels in the neutral medium with Hubble ex- 
pansion, its frequency redshifts. For a photon with initial fre- 
quency offset Av, the scattering optical depth to a distance r 
is 6 
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where «hi and H are the neutral hydrogen number density and 
Hubble constant at the redshift in consideration, (3 = Hr/c 
and t c = n m a c c/H. In what follows, we mainly adopt (3 as 
the distance variable, which proves to be convenient. At high 

6 We reduce the problem to a calculation for photons traveling in a medium 
with Hubble expansion velocity field with fixed Hubble constant at a given 
redshift. Strictly speaking, this approximation is only accurate if the distance 
traveled is much less than the Hubble radius. Such a condition is slightly bro- 
ken for the largest scales in the calculation, ~ 20% of the Hubble radius (see 
Fig.ff), which would introduce slight distortions in the scattering rate profiles 
on these scales. However, for our main purpose, the scales in consideration 
are generally much smaller and the above approximation is always valid. 
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redshifts, the conversion from f3 to comoving distance d is 
simply 



d =1176,3 
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The parameter t c can be regarded as a redshift variable for 
a given cosmology. Formally, it is the optical depth to the 
Hubble radius seen by a photon with frequency twice the line- 
center frequency. For Lya photons, we have 
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where the mass fraction of helium is taken to be one quarter. 

Equation (J8j can be inverted to find the distance traveled by 
a photon of initial frequency offset Av/vq for a given optical 
depth r, 
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For the solution to approach (3 = Av/vq, the one correspond- 
ing to a Dirac So cross-section, the frequency offset Av jvo 
needs to be large and/or the parameter r c needs to be small. In 
the regime that Av/vq <§C t c , (3 is proportional to {Av/vq) 2 , 
which implies that on sufficiently small scales the distance 
where most scatterings occur deviates from the expectation 
with So cross-section. The dependence on r in Equation ( fTTT i 
also means that, for a given initial frequency offset, the place 
where most scatterings occur is no longer at a single distance 
and instead it has a distribution. 

The spherical symmetry and Hubble velocity field allow 
a simple Monte Carlo calculation of the distribution of Lya 
photons from continuum between Lya and Ly/3 as a function 
of distance from the central source. The procedure is as fol- 
lows: 

i. Start a photon from the center (/?, = 0) with frequency 
offset Av/vq drawn according to the slope of the UV 
spectrum of the central source. 

ii. Draw a scattering optical depth t according to the ex- 
ponential distribution. 

iii. Find the distance f3 of traveling before scattering ac- 
cording to equation (fTTT i. 

iv. Determine the position of scattering, (3r = (3j + (3. That 

where cos# is drawn 



is, (3 f = y/ffW 1 
uniformly between 



2j3ij3 cos ( 
-1 and 1. 



v. Av/vq *— Av/vq- (3 (i.e., redshifted) and /?,■ <— (3f. 

vi. Repeat ii.-v. until Av/vq approaches the width of the 
thermal core (a few times 10~ 6 for situations here). 
Then start from i. again until the desired number of 
photons have been drawn. 

The cosmology and redshift are fully encoded in a single pa- 
rameter t c [eq. ( TTOb l . In the following discussions, we will 
assume t c = 0.0188 for Lya photons, which corresponds to 
z = 25. In step iv., we simply assume the direction after scat- 
tering is isotropic. We have tested that a more realistic angular 
distribution, such as a dipole distribution, has little effect on 
the resultant spatial distribution of Lya photons. 
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FIG. 1. — Radial probability distribution function (PDF) of "continuum" 
Lya photons around a central source at z = 25, which is proportional to r 2 
times the Lya scattering rate PcontM- The top and bottom panels are for 
Lya photons with initial frequency offset Av/vq of 10~ 3 and 10~ 2 , respec- 
tively. The dotted curve, dashed curve, rilled circles, and open circles are 
the distributions after the first, second, 10th and 100th scattering. The solid 
curve is the asymptotic distribution before core scatterings happen, which 
is the same profile of the scattering rate for photons with the given initial 
frequency offset. If a Dirac &d function is adopted for the scattering cross- 
section, the distribution would be a spike at the right edge of each panel, with 
Hr/c = Av/v a . 



In Figure Q] we show the spatial distribution of Lya 
photons, Wk. a (r, v), after k scatterings for two initial fre- 
quency offsets Av = v—v a , based on the above Monte Carlo 
procedure. The distribution has been normalized so that 
J Air ^Wu.aift ,v)df3 =1. In each panel, the dotted curve, 
dashed curve, filled circles, and solid circles are the distri- 
butions around the central source after the 1st, 2nd, 10th and 
100th scattering, respectively. As photons redshift towards the 
line center, the scatterings become more frequent and the dis- 
tribution approaches an asymptotic one, W a (r,v) (thick solid 
curve). Hence, in practice, to derive the total scattering rate 
of Lya photons, it is sufficient to perform the Monte Carlo 
simulation for the first few scatterings. 

If the Lya cross-section is assumed to be a Dirac So func- 
tion, for an initial Av/vq photon, all the scatterings would 
happen at a distance (3 = Av/vq. However, this is not the 
case if the frequency dependence of the cross-section is taken 
into account, as shown in Figure Q] (thick solid curves). The 
position where most of the scatterings happen is broadly dis- 
tributed at distances smaller than (3 = Av/i/q. Furthermore, 
the peak of the distribution depends on the initial frequency, 
the smaller the initial frequency offset AvJvq, the farther 
away of the peak from the position expected from a So func- 
tion cross-section. For example, as shown in Figure [T] for 
Av/i/q = 10~ 2 , the peak of the distribution is at a distance 
~ 45% of that expected from a So function cross-section, 
while for Av/vq = 10~ 3 it is only ~ 15%. As already men- 
tioned, this dependence on initial frequency can be under- 
stood by considering the first scattering with equation ( fTTT ): 
for a given optical depth r, (3 becomes increasingly smaller 
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than Av/vq as Av/vq becomes smaller. Equation ( fTTT ) also 
shows that for smaller t c , the scattering position would be- 
come close to that expected from a So function cross-section, 
which we discuss more about in the next section. 

To get the total scattering rate of the "continuum" Lya pho- 
tons for a general UV source, one needs to integrate over the 
frequency range between Lya and Lyp 



comoving distance (h ^pc) 
0.1 1 10 10 2 



^cont(r) = T GP S a / W a (r,i/)L v dv, 
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where tqp is the Gunn-Peterson optical depth to Lya photons, 
S a is the backreaction correction factor [eq. dj)], and L v is the 
luminosity of the source in terms of number of photons per 
unit frequency per unit time. 

The above dependence of the scattering position on the ini- 
tial frequency implies that the radial profile of the scattering 
rate would be steeper than the 1 jr 1 drop. With the Monte 
Carlo technique, we calculate the scattering rate profile from 
"continuum" photons, assuming a fiat UV spectrum of the 
central source ("fiat" here means equal number of photons 
per unit frequency per unit time). In Figure|2] the dashed blue 
curve is the result from the Monte Carlo simulation, while the 
dotted blue curve is the corresponding curve of the 1/r 2 drop 
with the same normalization of the UV spectrum of the central 
source. The true profile is steeper than 1/r 2 and the amplitude 
can be an order of magnitude higher at small distances. 

Close to the source the scattering rate scales as P CO nt(r) oc 
r -7/3 rp n j s behavior can be explained by the following ar- 
gument. The mean free path / of the photon emitted with 
frequency offset Av (i.e., the average distance between sub- 
sequent scatterings) scales as / oc (Ai/) 2 feq-dTTb in the limit 
of Af/Vo «t c ~ 10~ 2 ]. The change of frequency offset be- 
tween two subsequent scatterings is simply proportional to 
I. Therefore the number of scatterings that would occur be- 
fore the photon frequency moves closer to resonance and its 
mean free path drops significantly, scales as N sc oc Av/l oc 
(Ai/)" 1 . Since the photon direction changes almost at ran- 
dom after each scattering, the total distance it travels until it 
redshifts into the resonance scales as r oc ,/AQ oc (Avfl 2 . 
Conversely, we can say that photons reaching the resonance 
within distance r from the source, are emitted with frequency 
within Au oc r 2 / 3 from Lya resonance. Since the total num- 
ber of photons within this frequency range increases in di- 
rect proportion with Av, the photon scattering rate, which 
is proportional to the number of photons reaching the reso- 
nance within some region divided by its volume, scales as 
P CO nt(r)oc Ais/P (XT 7 ? 3 . 

3. RADIATIVE TRANSFER - "INJECTED" PHOTONS 

In addition to photons originally emitted between Lya and 
Ly/3 frequencies, which redshift into the Lya resonance, Lya 
photons are also produced by splitting of photons originally 
emitted between Ly7 and Lyman-limit frequencies. These 
photons will be first scattered as they redshift into the clos- 
est resonance. In case, following their absorption by a hy- 
drogen atom, the excited electron cascades directly to the 
ground state, the original photon would be re-emitted. Alter- 
natively, the electron can cascade via some intermediate state, 
in which case the original photon is split into several photons. 
Depending on the path of the cascade, the cascade products 
may include the Lya photon. The fraction of Lya photons 
made up by casca de of high resonances is typically less than 
~ 15% of the total dHiratal2006tlPritchard & F urlanett ol2005 
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FIG. 2. — The radial profile of Lya scattering rate around the central source 
at z = 25 as a sum of contributions from "continuum" and "injected" photons. 
In the top panel, the blue (dashed), magenta (thin solid), and black (thick 
solid) curves are those from "continuum" photons and "injected" photons 
and the total. The corresponding dotted curves are obtained assuming the 
scattering cross-sections for Lyman lines to be Dirac So functions. In the case 
of the "injected" photons, the dotted curve is almost on top of the magenta 
(thin solid) curve. Dashed vertical lines mark the maximum distance reached 
by Lyman series photons with the first five series labeled. In the bottom 
panel, the fractional differences in the scattering rates with respect to the 
case of So cross-sections are plotted for the "continuum" photons (dashed 
blue), "injected" photons (thin magenta, near 0), and the sum (thick black), 
respectively. A flat UV spectrum is assumed (see the text). 



IChuzhov & Sha piro 2007). However, as the distance traveled 
by photon before it redshifts into the closest np resonance is 
of order (v n +i/v n — l)c/H oc n~ 3 , the "injected" Lya photons 
are produced within much smaller volume and within ~ 10 
Mpc (comoving ) from the source o utnumber the "continuum" 
Lya photons dChuzhov et alj|2006l) . 

Since the "injected" Lya photons are injected directly into 
the resonance (hence their name), to derive their scattering 
rate profile it is sufficient to follow the path of their high res- 
onance progenitors. The Monte Carlo simulation is similar to 
what we perform for "continuum" photons, but we need to use 
cross-sections for higher resonance lines [i.e., different r c and 
Vd parameters in eq. ( TTTb for different lines] and record the 
probability of Lya production at eac h scattering. In the cal- 
culatio n, we use the data in Table 1 of |Pritchard & Fur lanettol 
(2006) for the probability of decay from the np state to the 
ground Is state and the probability of producing a Lya pho - 
ton through cascades from the np level (also see Hirata 2006). 
The scattering cross-section for higher resonance Lyman lines 
are calculated from oscillator streng th and E i nstein A coeffi- 
cient for np — > Is transition listed in M ortonl d2003l) (and ex- 
trapolations are used for higher n). 

The total scattering rate of Lya photons produced by cas- 
cade from np state with initial frequency v = v n + Av is 



Pmi,n(r) = T G pS a 
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where Wk, n (r, v) is the normalized spatial distribution for Lyn 
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FIG. 3. — Radial distribution of "injected" Lya photons around a central 
source at z = 25, which is proportional to r 2 times the Lya scattering rate 
Pi n j(r). Shown here are "injected" Lya photons from Ly<5 (top panel) and Lye 
(bottom panel) resonances with initial frequency offset 10""*. In each panel, 
the dot-dashed, long dashed, short dashed, and dotted curves are contributions 
from Lyn (n = S,e) photons that experience no more than 1, 2, 10, and 20 
scatterings before destroyed to produce Lya photons. The thick solid curve 
is the total contribution from all Ly« photons. The distribution from a Dirac 
So cross-section would be a spike at the right edge of each panel, with Hr jc = 

ICr 4 . 

photons that experience k scatterings, /d es ,n is the probability 
that the original Lyn photon is destroyed, f a n is the chance 
that the destruction of the Lyn photon results in production 
of the Lya photon (i.e., that the electron cascade goes via 2p 
rather than 2s state), and other symbols have the same mean- 
ing as in equation ( TTSl i. Since at each scattering there is a 
significant probability, 0.12 < /d es ,n ^ 0.23, that the original 
Lyn photon is destroyed, at large values of k, (1 -/des.n)* 1 ap- 
proaches zero. That is, the total distribution of "injected" Lya 
photons can converge by considering those np photons that 
experience only a small number of scatterings. 

Figure [3] illustrates the evolution of Lya photon distribu- 
tion for higher resonance photons (Ly<5 and Lye, with initial 
frequency offset 10~ 4 ) that experience no more than k scatter- 
ings (£=1,2, 10, and 20) before destroyed. As expected, the 
distribution approaches quickly to the total distribution (solid 
curve), similar to the "continuum" case. At the same initial 
frequency offset, higher resonance photons lead to a peak po- 
sition closer to that expected from <5d function cross-section. 
The reason is that a higher resonance line has a smaller cross- 
section [thus a smaller r e parameter in eq. dTTbl . 

In contrast to "continuum" photons that typically reach the 
resonance after multiple scatterings, a photon emitted with 
frequency blueward of a higher Lyman resonance can produce 
a resonant Lya photon just after a single scattering. If such 
a photon has mean free path I oc (A^) 2 as the "continuum" 
photon case (see § 2), the distance it travels before produc- 
ing a Lya photon is just r oc / oc (Az/) 2 and one would expect 
the scattering rate profile Pmj. n (r) oc Av/r 3 oc r~ 5/2 . However, 
this can only formally happen extremely close to the central 
source [Au/u„ <C r c , see eq.dTTbl because of the fast drop of 



the corresponding r c , and proximity effect would take over at 
such scales. Since both the oscillator strength / and the spon- 
taneous decay rate 7 drops as n~ 3 at large n, the characteristic 
scattering optical depth t c in equation ( fTTT i drops as n~ 6 . Even 
for n = 3, t c has dropped to a 10~ 5 level. The low value of 
t c means that, at a given initial frequency offset, the scatter- 
ing rate distribution for injected photons is almost the same 
as that expected from 5q function cross-section with the peak 
position approaching (3 = Av/v n . For scales of interest (i.e., in 
the regime of Av jv n > r c ), the mean free path / oc Av instead 
of {Av) 2 for higher resonance photons, thus each individual 
P'mi.n(r) follows the 1/r 2 profile. 

The solid magenta curve in Figure [2] shows the Lya scat- 
tering rate profile from injected photons for the same flat UV 
spectrum of the central source as in the "continuum" case. 
Each of the steps in the curve reflects the place that a higher 
resonance line starts to contribute to produce injected Lya 
photons. At a given initial frequency offset, the scattering 
rate of injected Lya photons is dominated by those produced 
by higher level resonance lines. As mentioned above, the lo- 
cation of the injected Lya photons from these resonance lines 
is close to the expectation of <5o-form cross-section as the cor- 
responding t c decreases fast. Therefore, the true scattering 
rate profile is quite close to that from fe-form cross-section, 
which is a sum of a series of truncated 1/r 2 functions (dotted 
magenta curve in Figure [2] almost on top of the solid magenta 
curve). At a distance j3 = 10~ 5 , the true profile is only higher 
by ^4% in this particular case. 

Although each individual P m i„{r) follows a truncated 1 /r 2 
profile, the overall profile from summing over all "injected" 
photon contributions no longer follows a 1/r 2 profile because 
the truncation place for the individual component depends on 
n. It can be shown that at small scales, the overall profile 

Y Jn p H,ni r ) °c r ~ 7/3 dChuzhov et al.1 120061) . interestingly the 
same dependence as the "continuum" photon case (see § 2). 
To see this, we note that, with r oc (v- v„)/v m the individual 
component P m i n {r) is simply oc L u dv /(4irr 2 dr) oc v n fr 2 for 
a flat source spectrum L v = constant (in photon number per 
unit frequency per unit time). The slope of the overall profile 
is then -2 + dln(J2„ v n )/d\nr. Since v n oc 1 — 1 /n 2 , we have 
J2„ v n °c n and r oc 1 /n 3 in the large n limit. Therefore, the 
overall profile ^2 n Pmj. n (r) oc r" 7 / 3 , which is what is seen in 
Figure |2 

4. THE TOTAL SCATTERING RATE 

The total scattering rate of Lya photons is obtained after 
summing over all Lyman series: 

00 

PqSX) ~ ^cont 

(r) + ^> inj ,„(r). (14) 

The solid black curve in Figure [2] shows the profile of the to- 
tal scattering rate. For the case shown here (z = 25), "contin- 
uum" ("injected") photons dominate at distances to the cen- 
ter greater (less) than /?~2x 10~ 3 (~ 3Mpc comoving). At 
smaller distance, even though the amplitude from "contin- 
uum" photons can be an order of magnitude higher than the 
naive rate given by the fe-form cross-section, the domina- 
tion of "injected" photons reduces the difference in the total 
scattering rate. As for the slope, the overall scattering rate 
is slightly steeper than the sum of a series of truncated 1/r 2 
drops (dotted black curve). To the first order, the overall scat- 
tering rate is about 40% higher than the naive calculation in a 
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wide range of distances (see the thick black curve in the lower 
panel). 

The calculation of the scattering rate profile becomes much 
easier by assuming a <5o-form scattering cross-section. From 
the above example, we see that to the first order, the exact 
profile that takes into account the frequency dependence of 
the cross-section and the radiative transfer effect can be ob- 
tained by applying a correction to that from the calculation 
with So cross-section. This correction is almost a constant, 
increasing the amplitude by a few tens of percent. Apparently, 
the correction factor depends on redshift, larger at higher red- 
shift because of an increase in the r c parameter [eq. (ITTT ll. We 
perform simulations at a series of redshifts and find that the 
fractional difference f c between the exact and the naive calcu- 
lations of the scattering rate can be approximated as a constant 
at a given redshift, 



comoving distance (h ^pc) 



0.1 



f c = 0.381 



l+z 
26 



0.66 



(15) 



This approximation underestimates (overestimates) the cor- 
rection at small (large) scales. For example, at z = 25, / e =44% 
(32%) at (3 = 10~ 5 (10~ 2 ). A more accurate fit that accounts for 
such a tilt is 

/c = 0.364 (i±£) a71 + (0.0115^-0.0443)log^3, 

(16) 

which works almost perfectly for 15 < z < 35 and 10 5 < (3 < 
10~ 2 (roughly corresponding to 0.01-20 /r'Mpc comoving). 
The first term on the right hand side always dominates and the 
second term accounts for the slight difference in the slope of 
the two profiles. With this fitting formula, the profile from the 
naive calculation can be corrected by multiplying the scatter- 
ing rate by a factor of 1 +f c . 

All the above examples assume a flat UV spectrum, L v oc 
v a ~ x with a s = 1 between Lya and Lyman limit, where L v is 
the luminosity in terms of number of photons per unit fre- 
quency per unit time. The spectrum slope a s depends on 
the nature of the central source. For example, for Popula- 
tion III stars, a s =1 .29, while for Population II stars a s = 0. 14 
(Barkana & Loeb 2005). Given the narrow frequency range 
between Lya and Lyman limit, we do not expect a strong de- 
pendence of the scattering rate profile on a s . In the case of 
adopting the So cross-section, the scattering rate for an in- 
dividual component is simply oc L u dv /{An^dfJ) oc Vq '(1 + 
(3) a ~ l I ' 1 , where the relation (3 = (v-vq)/vq is used. The de- 
parture from the 1 jr 1 drop with respect to the fiat spectrum 
is introduced through and (1 + (3) a ~ l and both are small 
given the narrow frequency range of Lyman series and (3 <C 1 . 

We perform simulations for two cases with quite different 
spectral slopes, a s = 2 and a s = 0. In Figure |4] we plot ratios 
of scattering rates for the two cases. The central source lu- 
minosity is normalized so that the numbers of photons emit- 
ted between Lya and Lyman limit for the two cases are the 
same. That's why the a s = 2 case has a higher rate from "in- 
jected" photons (thin solid curve) but lower rate from "con- 
tinuum" photons (dashed curve). The ratios for the individual 
components are similar to the 5q cross-section results (dotted 
curves). However, the steepening of the slope of the scatter- 
ing rate for the "continuum" photons at small scales, which 
increases its contribution to the total scattering rate, reduces 
the difference in the total scattering rates for a s = 2 and a s = 
(thick solid curve). Even with the large difference in the spec- 
tral slope, the difference in the overall scattering rates is only 
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FIG. 4. — Fractional differences in the scattering rate profiles for different 
slopes of central source spectra at z = 25. The luminosity of the central source 
(number of photons per unit frequency per unit time) between Lya and Ly- 
man limit is assumed to follow L„ oc u"'^' with a s = 2 and 0, respectively. 
The spectra are normalized to have the same number of photons between Ly a 
and Lyman limit. The dashed and thin solid curves are for rates from "con- 
tinuum" and "injected" photons, respectively, and the thick solid curves is for 
the total rate. Dotted curves are for the corresponding cases with Dirac 8d 
cross-sections. 

at a level less than 10% for (3 < 10~ 2 . A comparison between 
cases using spectral slopes of Population II (a s = 0.14) and 
Population III (a., = 1 .29) stars shows a pattern similar to what 
is seen in Figure|4]with the difference in the overall scattering 
rates being at a level of 5%. 

5. SUMMARY AND DISCUSSION 

We investigate the effects of frequency dependence of scat- 
tering cross-section and radiative transfer on the scattering 
rate profile of Lya photons around a central UV source, which 
is relevant for the pumping of the 21cm line in the high red- 
shift universe. 

Because of the frequency dependence of the scattering 
cross-section, a photon between Lya and Ly(3 frequency 
("continuum" Lya photon) experiences a small number of 
wing scatterings before it redshifts to the core frequency and 
starts core scattering. Core scatterings, which determine the 
scattering rate, happen at distances much less than that ex- 
pected from a So scattering cross-section. For scatterings of 
higher resonance photons that produce "injected" Lya pho- 
tons, the effect is weak owing to the fast drop in the cross- 
section. We have shown that for a single UV source the 
scattering rate profile of "continuum" Lya photons is signif- 
icantly steeper than l/r 2 (previously expected from 5d scat- 
tering cross-section), approaching r~ 7/3 near the source. The 
scattering rate profile of "injected" Lya photons from an indi- 
vidual high Lyman resonance closely follows a truncated l/r 2 
profile, while the overall profile from all higher Lyman reso- 
nances coincidentally also approaches r~ 7 ' 3 near the source. 
At 15 < z < 35, the total scattering rate from "continuum" 
and "injected" photons is 30-50% higher than the naive cal- 
culation that adopts Sq cross-sections, on scales of ^0.01- 
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20/i _1 Mpc (comoving). We also find that, when radiative 
transfer effects are properly accounted for, the slope of the 
scattering rate profile does not have a strong dependence on 
the spectral slope of the central source. 

In our calculations we have assumed that during the early 
stages of reionization the gas temperature changes adiabati- 
cally. If instead, the gas was significantly heated, the Doppler 
core for Lya photons would be increased and their scattering 
rate profile would be even steeper. However, since the temper- 
ature of neutral hydrogen does not exceed ~ 10 4 K, the scales 
above ~ 100 kpc (comoving) would be virtually unaffected. 

As a consequence of the higher scattering rate with respect 
to the naive calculation, the spatial fluctuations in the pump- 
ing radiation field produced by multiple UV sources would 
be significantly higher as well. During the early stages of 
reionization, when the UV intensity is relatively low, the dif- 
ferential brightness of the 21cm signal scales almost linearly 
with the scattering rate of the Lya photons and the corrections 
to the scattering rate translate into similar corrections to the 
21cm brightness. Such a correction would increase the size 
of the expected Lya spheres around first sources (e.g., ICenl 



120061: IChen & Miralda-Escude1l2006l) . Moreover, if the radia- 
tive transfer effect is overlooked, the fluctuation power spec- 
trum estimated from 21cm observations would be skewed on 
scales up to a few tens of Mpc, which would lead to inaccu- 
rate inference of the matter fluctuation power spectrum (such 
as the amplitude and the running of the spectral index). There- 
fore, until the epoch when the intensity of pumping radiation 
reaches a saturation level (if such epoch in fact exists), the 
correct interpretation of the 21cm signal requires taking into 
account radiative transfer (mainly for "continuum" Lya pho- 
tons), as described in this paper. 
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